Heating mechanisms in radio frequency driven ultracold plasmas 
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Several mechanisms by which an external electromagnetic field influences the temperature of a 
plasma are studied analytically and specialized to the system of an ultracold plasma (UCP) driven 
by a uniform radio frequency (RF) field. Heating through coUisional absorption is reviewed and 
applied to UCPs. Furthermore, it is shown that the RF field modifies the three body recombination 
process by tunnel-ionizing electrons from intermediate high-lying Rydberg states, whereby the three 
body recombination rate and associated plasma heating is suppressed. Heating through coUisionless 
absorption associated with the finite plasma size is calculated in detail, revealing a temperature 
threshold below which coUisionless absorption is ineffective. It is found that, for large RF amplitude, 
the ponderomotive potential replaces the electron temperature in functional dependencies of several 
plasma quantities. 

PACS numbers: 52.55.Dy, 52.27.Gr, 52.50.Qt, 34.80.Lx, 52.50.Sw, 36.40.Gk 
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I. INTRODUCTION 



Conventional plasmas are formed when atoms are ion- 
ized by strong electric fields or collisions with other par- 
ticles. Due to the large excess energy inherent in such 
ionization processes, the resulting electron temperature is 
typically comparable to the ionization potential, which is 
on the order of an electronvolt, equivalent to some 10'* K. 
In marked contrast, ultracold neutral plasmas (UCPs), 
created by photo-ionization of a cloud of laser-cooled 
atoms have an electron temperature close to 1 K. 
UCPs typically consist of some 10^ singly-ionized atoms 
localized in a millimeter-sized cloud of Gaussian density 
profile, with a correspondingly low particle density 2]. 
The combination of low temperature and low density 
makes UCPs unique plasma systems. They can be close 
to the strongly-coupled regime where the Coulomb inter- 
action energy between the particles exceeds the thermal 
energy, as is quantified by the coupling parameter 



diagnostics. In contrast, UCPs evolve on the time 
scale of picoseconds to microseconds. This enables 
excellent time-resolved diagnostic techniques, including 



F = 



47reo rujfcsT 



(1) 



exceeding unity, where e is the electron charge, eo 
the vacuum permittivity, fcs Boltzmann's constant, T 
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the plasma temperature, and r^, — [3/(47rn)] ' with 
n the number density. Due to their high coupling 
parameter, UCPs behave in many respects similar to 
strongly-coupled plasmas near solid state density, such 
as laser- ionized atomic clusters Q or thin films 
inertial confinement fusion targets Q and astrophys- 
ical plasmas @. The dynamics of solid state density 
plasmas, however, takes place at the time scale of the 
inverse plasma frequency, which lies in the attosecond 
to femtosecond regime. This seriously complicates 



charge particle detection 
fiuorescence monitoring 
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, absorption imaging [8| and 
j. In addition, the careful 
preparation and ionization of atomic clouds allows 
accurate control over the initial temperature, density 
profile, and ionization state. UCPs may therefore serve 
as versatile and experimentally accessible model systems 
for high-density plasmas that are difficult to diagnose. 

An important class of experiments on solid state 
density plasmas involves plasmas created by laser 
irradiation of atomic clusters in a gas jet. Characteristic 
of these experiments is that the laser pulse length is 
comparable to the lifetime of the plasma. Therefore the 
studied system typically consists of a cluster plasma that 
is not only near to strongly-coupled, but is also strongly 
driven by a radiation field. This leads to complicated 
dynamics that is difficult to unravel Q. Research on 
laser-cluster interaction would therefore benefit from 
UCP experiments in which this interaction is mimicked. 
Since atomic clusters are typically smaller than the laser 
wavelength, the appropriate model system is an UCP 
driven by a strong radio-frequency (RF) field. Interpre- 
tation of observations in such experiments on RF driven 
UCPs, however, requires a detailed understanding of the 
mechanisms by which the RF field and the UCP interact. 
In this paper, we consider how the RF field influences 
the plasma temperature, both directly through RF 
energy absorption mechanisms and indirectly through 
modification of the three body recombination process, 
the latter being a main heat source in UCPs. 

In current UCP experiments, RF fields are used in a 
diagnostic way to probe plasma modes. Plasma reso- 
nance can be detected as an increased yield of electrons 
leaving the UCP llOll . Combined with knowledge of the 
mode properties |11| . this can be used to determine the 
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plasma density and expansion as a function of time. 
Using the same technique, the presence of acoustic or 
Tonks-Dattner modes in an UCP has been observed 
in addition to the fundamental mode [l^l- In these 
experiments, the collective response of the plasma 
electrons to the RF field has been studied in quite 
some detail. However, the RF amplitude is kept low to 
avoid disturbances other than plasma resonances, and 
little attention is paid to other interaction mechanisms. 
Nevertheless, as we will describe in this paper, the RF 
field influences the plasma also via incoherent processes. 
In their Tonks-Dattner modes experiment, Fletcher 
et al. dl] indeed observe the onset of field-induced 
effects at large probing amplitudes. Although lower RF 
amplitudes justify the use of standard plasma quantities, 
such as the Spitzer collision frequency applied in the 
interpretation of the fundamental plasma resonance 
measurements [11], or the Debye length mentioned in 
support of the analysis of the Tonks-Dattner modes 
[12| . one should be aware of the possible high-amplitude 
modifications of such quantities induced by the RF 
field. Finally, the expansion of an UCP is driven by 
the thermal pressure of the electrons. It is there- 
fore important to understand the various ways in which 
the RF field contributes to the heat budget of the plasma. 



This paper is organized as follows. We consider two 
mechanisms by which the UCP can directly absorb en- 
ergy from the RF field: coUisional absorption and colli- 
sionless absorption due to the finite size of the plasma. 
The first of these has been studied extensively already in 
other contexts [T3l - [23j . In Section we therefore only 
cite the main results from literature and discuss their 
relevance for RF driven UCPs. Roughly speaking, for 
high RF amplitude it is found that in several kinetic 
plasma quantities the thermal velocity is effectively re- 
placed by the quiver velocity induced by the RF field, 
or equivalently that the temperature is replaced by the 
ponderomotive potential. This suggests more generally 
that other processes that depend on the temperature are 
analogously modified by the RF field. In Section IIII[ 
we confirm this notion for the process of three body re- 
combination, and show that the recombination rate, and 
associated plasma heating, is suppressed by the RF field. 
Next, in Section ITVl we consider the coUisionless absorp- 
tion mechanism mentioned above, which has been mainly 
studied in the context of solid-state density plasmas [2J- 
Isij . We show that the approximations usually adopted 
are not appropriate for UCPs. We provide an improved 
description by specializing a derivation of the coUisionless 
absorption rate due to Zaretsky et al. [2^ to the case of 
UCPs. We conclude and summarize in Section IVl 



II. COLLISIONAL ABSORPTION 

A. Collision frequency 

At low to moderate RF field strengths, the energy ab- 
sorption of a plasma is dominated by collisional absorp- 
tion, or inverse Bremsstrahlung [3^ . The physical cause 
of the absorption is that individual electrons, oscillat- 
ing due to the RF field, defiect in the Coulomb fields of 
the approximately stationary ions, resulting in a net en- 
ergy gain. The average effect of the Coulomb fields can 
be described phenomenologically as an effective frictional 
force F = —mVciV in the equation of motion of the elec- 
tron, and the energy absorption rate per electron by the 
power Pei = —{F-v). Here, m is the electron mass, 
is the effective electron-ion collision frequency, and v is 
the electron velocity. Expressing the velocity in terms of 
the driving electric field gives |13{ 



2iyciUp, 



(2) 



where Up = (ei?o)^/(4ma;^) is the quiver energy, or pon- 
deromotive potential, in the RF field with amplitude Eq 
and frequency uj. Here and in the remainder, we assume 
a linearly polarized RF field, and absorb any field en- 
hancement due to plasma resonance in the magnitude 
Eq. Importantly, Eq. ^ defines the collision frequency 
as merely a scaled absorption rate, rather than predicting 
the absorption from a predetermined collision frequency. 
Consequently, Uei is not necessarily equal to the Spitzer 
collision frequency [33j 



.r^/^iuA, 



(3) 



which is commonly used for plasmas without RF fields. 
Nevertheless, the collision frequency Eq. ^ is some- 
times used for driven plasmas as well, and also in the 
context of RF absorption by UCPs [11 |3||. In Eq. 

singly ionized atoms are assumed, LOp is the plasma 
frequency, and In A is the Coulomb logarithm that will 
be discussed below. 



Underlying any calculation of the collisional absorption 
rate is some model for the scattering of an electron by 
the Coulomb field of an ion, which generally depends on 
the electron velocity. Because two velocity scales are in- 
volved, namely the thermal velocity Uth = -y/ZcsTe/m and 
the quiver velocity magnitude Uosc — ^^^/{muj)^ the col- 
lision frequency depends on the ratio fosc/t'th- Here, Te 
is the electron temperature of the plasma. The effective 
collision frequency has been calculated first by classical 
kinetic theory using the Landau collision integral [l3l.[l3|. 
The result can be written as J 9.1 
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FIG. 1: Effective electron-ion collision frequency for colli- 
sional absorption scaled to the Spitzer collision frequency, as 
a function of the ratio of quiver velocity to thermal velocity. 
Solid black line: collision frequency Eq. (Q; dash-dotted blue 
line: weak field limit Vosc <C Vth given by Eq. ((5]); dashed red 
line: strong field limit Vosc S> Vth given by Eq. (O. 



B. Coulomb logarithm 

A second important consequence of the RF field is that 
the Coulomb logarithm In A in Eq. (jS]) is modified. This 
is particularly relevant for UCPs because the traditional 
expression In A = In (F"^/'^) looses its validity in case of 
strong coupling F > 1. The Coulomb logarithm arises 
from cutting off the Coulomb collision integral at both 
large and small impact parameters in elementary calcu- 
lations of the scattering cross section of an electron by an 
ion [36]. However, the physical arguments used to choose 
these cut-offs are traditionally based on thermal electron 
velocities only, and the cut-ofFs will change when in ad- 
dition the quiver velocity is taken into account. This 
can be confirmed by explicit calculation [17], yielding 

In A « In (femax/^min), with 



VcS 



max (w, ujp) ' 

„2 



(6) 

(7) 
(8) 



where_oF2 denotes the generalized hypergeometric func- 
tion \3&\ that has the limiting forms 



2F2(...) 




1.0 



Vosc < Vth 
Vosc > Vth- 

(5) 



More advanced and alternative calculations largely 
confirm these results Il5l420l. 



Here the classical limit < e^/(2eo^) has been as- 
sumed, where 2TTh is Planck's constant. Eqs. ([5][5]) show 
that also in the Coulomb logarithm, as before, the quiver 
velocity effectively takes over the role of the thermal ve- 
locity in the limit Uosc ^ Vth- This suggests more gener- 
ally that kinetic processes in UCPs that depend on the 
electron temperature may be strongly modified by the 
presence of an RF field. In the next section, we further 
validate this notion by showing that the three body re- 
combination rate in an UCP can be strongly suppressed 
by application of an RF field. 



The collision frequency of Eq. is plotted in 

Fig. [T] as a function of the velocity ratio. In RF 
experiments with UCPs, this ratio can vary over the 
full range t;osc ^ vth to Vosc ^ Vth [131 • The decrease 
of the collision frequency for increasing Wosc can be 
understood physically from the well-known fact that the 
Rutherford scattering cross section for an electron by 
an ion decreases proportional to w^"', so that driving 
the plasma stronger makes the electrons less susceptible 
to defiections and hence to energy gain. Note that the 
Spitzer frequency Eq. ^ with Eq. ([T]) substituted is 
proportional to v^^, while the second line of Eq. ([5]) 
contains the factor (I'th/wosc)^- Effectively, therefore, 
and apart from a logarithmic factor, the content of 
Eq. (U]) is that the thermal velocity is replaced by 
the quiver velocity in the collision frequency when 
Vosc ^ Vth- In fact, this effect is such that the col- 
lision absorption rate Pd given by Eq. ([2]) decreases 
with field strength as E^^ rather than increases, which 
is a well-known phenomenon in laser-plasma physics [3^ . 



III. THREE-BODY RECOMBINATION 

In the process of three body recombination (TBR), an 
electron recombines with an ion, while the excess poten- 
tial energy is carried away by a second electron. TBR is 
particularly significant in UCPs because the TBR rate 

R is proportional to Te according to conventional 
theory [33, HI]. Numerical simulations based on the 
Te dependence are indeed in agreement with UCP 
expansion measurements [3^ . As we will show now, 
the presence of an RF field, and the correspondingly 
increased velocity of the electrons, suppresses this 
divergent TBR rate. 

We determine the TBR rate along the lines of the 
calculation by Hinnov and Hirschberg [s^], adapted to 
the situation in which Vosc 3> Vth- Let us first briefly re- 
view the conventional case where the RF field is absent. 
Quantum mechanically, a TBR event may be described 
as an electron making a cascade of transitions between 
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adjacent energy levels of an atom until it reaches the 
deeply bound states. Under conditions applicable to 
UCPs, these transitions are mainly effected by collisions 
with other, free electrons. The process is illustrated 
in the left panel of Fig. [2j Considering an electron 
at any particular energy level Uk < 0, there is both a 
finite probability that the next collision will result in 
an upward transition, and a finite probability that a 
downward transition results. It can be shown [s^l that 
the upward transition probability increases with respect 
to the downward transition probability as grows 
closer to the continuum, and that upward transitions 
dominate for levels less than an energy fc^Te below 
the continuum. Any electron ending up in the energy 
band —ksTf. < U < 0, shown in gray in Fig. [51 is 
therefore likely to re-ionize, while electrons below this 
band are likely to fully recombine. Hence one may 
approximate the amount of recombining electrons with 
those electrons that skip the band altogether by making 
a direct collisional transition from the continuum to 
anywhere below the bottleneck level U* = —ksTe- The 
validity of this approximation has been confirmed by 
extensive Monte Carlo simulations [11]. Summing the 
probabilities of such transitions over all possible initial 
and final energies of the recombining electron and over 
all possible energies of the free electron, one finds indeed 
the usual TBR rate proportional Tg W^ - 



(d) (c) 
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FIG. 2: Energy diagram of three body recombination with 
and without RF field. The RF field induces a Stark shift, 
which is approximately equal to Up for the high energy levels 
(a) and much smaller for deeply bound states (b). One high 
energy level f/fc < is drawn. The gray bands show energies 
above the bottleneck level from which re-ionization is likely, 
(c): electron that re- ionizes after a collisional transition from 
the continuum to above the bottleneck level; (d): electron 
that recombines by making a cascade to deeply bound states 
after a collisional transition to below the bottleneck level. On 
the sides the energy scales as used in the text are indicated. 



When an RF field is present, essentially two things 
change in this picture, as illustrated by the right 
panel of Fig. O First, the energy of both free and 
bound electrons will change due to the field. For 
free electrons, the energy increment is just the quiver 
energy Up = mvl^^/A. For bound electrons, the energy 
increment is the AC Stark shift. However, the energy 
levels that are significant for the TBR rate are the 
high-lying Rydberg states, for which the AC Stark 
shift is approximately equal to Up as well |40|. This 
corresponds to classical Kepler orbits with approximate 
free electron quiver motion superposed. This uniform 
Stark shift is shown in Fig. Neglecting the low-level 
deviations from the uniform Stark shift, and placing 
the origin [/ = at the shifted continuum boundary, 
the situation would be energetically equivalent to the 
case without RF field if that was the only RF effect. 
Secondly, however, the RF field ionizes atoms by the 
process of tunnel ionization. The probability of such 
ionization from energy level Uk is exponentially small 
if Uk is more than an energy ~ Up below the (Stark 
shifted) continuum, but is very large if Uk is higher 
than that [4l|. This, completely analogous to the case 
without RF field, defines an energy band just below 
the continuum from which electrons are much more 
likely to be re-ionized than to proceed with a downward 
collisional cascade, although now the energy band has 
widened to —Up < U < 0. This is shown as the thickened 
gray band in Fig. [51 Thus the primary infiuence of 
the RF field on the TBR process is that the bottleneck 



level is lowered from U* = —ksTe to U* = —Up, replac- 
ing once more the thermal velocity by the quiver velocity. 



Taking account of the adapted bottleneck level, the 
TBR rate can be evaluated as done by Hinnov and 
Hirschberg. By the method of detailed balance, under the 
hypothetical condition of thermal equilibrium the rate of 
collisional transitions from the continuum to the energy 
level Uk < is equal to the rate of the inverse process, 
which is ionization from level Uk caused by electron im- 
pact. From the well-known |37| cross section Sk{U) for 
a collisional energy transfer of at least \Uk\ from a free 
electron with energy f7 to a bound electron with energy 
Uk, the rate of ionization from level Uk per unit plasma 
volume is 

poo 

Rk^ / nknevSk{U)f{U)dU. (9) 

Here, f{U) is the energy distribution function of free elec- 
trons, V — \/lU jm, and rife is the density of electrons in 
level Uk- The rate of TBR via level fc, which is the in- 
verse process, is obtained by substituting for rife the equi- 
librium value from the Saha equation [SQl, because the 
two rates must be equal at equilibrium. Evaluating the 
integral in Eq. ^ assuming a Boltzmann distribution 
f{U), substituting the Saha value for nk, and summing 
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over all energy levels below U* , gives the total TBR rate 



R — ^ Rk 



U' 



F 



Uk \ gkdUk 



^ 8g.e§m2 (fc^Te)V-oo \kBTeJ dUk/dk' 

(10) 

where F{x) = exp(— x) Ei(a;) — 1/x with Ei the exponen- 
tial integral [ill, gk is the degeneracy of level k, and gi 
is the degeneracy of the ionic ground state. The sum 
over states has been approximated by an integral over fc, 
after which the integration variable has been changed to 
the energy Uk- So far, the RF field has not entered the 
discussion, so that Eq. ([TU|) holds whether an RF field is 
present or not. In the case with field, the field effects con- 
veniently separate, that is, the bottleneck modification is 
expressed in the integration limit, while the Stark shift 
structure is fully contained in the factor gk/{dUk/dk). 
Noting that \Uk\/{kBTe) ^ 1 over the whole integra- 
tion domain, the function F can be approximated by its 
asymptotic value F ^ {ksTe/Uk) ■ Then the remaining 
integral does contain the field effects, but is independent 
of the temperature, so that the temperature scaling of the 
TBR rate derived here is insensitive to our simplifying 
assumptions of a sharply defined tunnel ionization band 
and a uniform Stark shift. Assuming Rydberg atoms, for 
which Uk — ~Ry/k^ with Ry — 13.6 eV the Rydberg 
energy, and gk ~ k^, integration of Eq. pUj) gives 
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assuming g^ = 2. Within a factor of order unity, this 
three body recombination rate is equal to the accepted 
result for the case without RF field [3^, except that 7/2 
powers of fc^Te have been replaced by the quiver energy, 
or equivalently that 7 powers of the thermal velocity 
have been replaced by the quiver velocity. This reduces 

9/2 

the strongly divergent behavior R (x Te to the much 
milder dependency R oc T^^. Thus three body recom- 
bination, and the associated plasma heating, may be 
significantly suppressed by the application of an RF field. 



IV. COLLISIONLESS ABSORPTION 
A. Absorption models 

Even without the presence of electron-ion collisions, 
individual electrons in a plasma can absorb energy 
from an applied electric field. For bulk plasmas, this 
coUisionless absorption effect is the well-known Landau 
damping [3^, in which electrons can gain net energy 
from a high-frequency propagating electric wave, despite 



the fact that the high-frequency electric force tends to 
cancel out on the average. This is possible when the 
thermal velocity of the electron is close to the velocity 
of the wave, so that the electric field is approximately 
static in the electron frame of reference. Essential for 
this mechanism is a correlation between thermal motion 
and applied field. In plasmas of finite size, such as an 
UCP, the thermal motion of electrons is necessarily 
confined by the plasma boundaries, so the assumption 
of rectilinear motion implicit in the Landau damping 
mechanism of bulk plasmas is no longer appropriate. 
Rather, the electrons perform quasi-periodic motion in 
the electrostatic potential of the plasma, as is detailed 
below. Furthermore, the electric field in the plasma 
is homogeneous rather than a propagating wave when 
the applied wavelength is much larger than the plasma 
size, such as in the case of an RF field applied to an 
UCP. Nevertheless, electrons may on the average gain 
energy, and this is again due to a correlation between 
the thermal motion and the applied field. This is why 
the coUisionless absorption of finite plasmas has been 
called Landau damping as well [1^ |4^, although the 
character of the correlation is quite different. In this 
section, we calculate the RF energy absorption of an 
UCP by this mechanism. 

First, we mention a number of other approaches to 
coUisionless absorption and argue why these are less ap- 
propriate for UCPs in RF fields. In the above descrip- 
tion of coUisionless absorption, the applied field plays 
the role of a perturbation on the thermal motion of the 
electrons. One may change perspective and look at the 
quiver motion of the electrons as being the primary mo- 
tion, perturbed by a thermal one. Because the details of 
the thermal motion are determined by the details of the 
plasma potential, this can be interpreted as an oscillat- 
ing electron having interaction with the plasma potential 
itself. This view is particularly appropriate when the po- 
tential can be approximated by an infinitely deep well, so 
that the 'interaction with the potential' simply becomes 
'collisions with the plasma boundary'. Then the collision 
frequency of electrons with the plasma boundary is on 
average 



V 

a 



(hard wall model). 



(12) 



where a is the plasma size, and v is the characteristic 
velocity of the electrons that is taken to be the thermal 
velocity (28j , a combination of thermal and quiver veloc- 
ity [23] or Fermi velocity [25)] depending on the model 
used. On average the electrons gain an energy 2C/p per 
hard wall collision, in analogy with Eq. The re- 

sult also follows as a special case from the more 
general Landau daniping approach when specialized to a 
hard wall potential [28| . While a flat potential with hard 
walls, and hence the resulting absorption rate 2vC7p/cr, 
may be a good approximation for large metallic clusters 
|25l [27} ■ it is not for UCPs. In the process of creation 
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of an UCP from an atomic cloud, part of the electrons 
escape from the plasma immediately after photoioniza- 
tion of the cloud. This continues until the accumulated 
charge imbalance self-limits further loss of electrons. The 
resulting spherically symmetric Coulomb potential of the 
UCP with a typical Gaussian density distribution is Q 



U{r) = C/o 



1 



2r 



■ erf 



(13) 



where erf(r/cr) denotes the error function [42], and r 
is the distance to the cloud center. The depth of the 
potential saturates to Uq ^ ksT^. by nature of the 
charging process. Clearly, the hard wall potential is 
not a very good approximation in this case and a more 
detailed calculation of the energy absorption is necessary 
to account for the smoothness of the potential. 

Another absorption mechanism that is considered 
important for large metal clusters is the Brunei effect 
[26| . in which electrons at the plasma boundary are 
pulled out of the plasma by the applied electric field and 
then driven back into the plasma as the field reverses 
direction. When the plasma is sufficiently overdense, the 
interaction effectively stops once the electron has moved 
deeper into the plasma than the skin depth, resulting in 
net energy gain because the electron cannot be brought 
back to rest by the evanescent field. The resulting 
absorption rate, divided by 2Up for comparison, gives 
again the hard wall collision frequency Eq. (fT2)) . with 
V the high-frequency velocity. In an UCP, however, the 
Brunei mechanism is not in effect either, since typically 
the skin depth, which is comparable to c/ujp with Up the 
plasma frequency, is much larger than the plasma size. 

Finally, when the applied field is so strong that the 
oscillation amplitude of individual electrons is compa- 
rable to or larger than the plasma size, one can hardly 
speak of the applied field as a perturbation, and other 
descriptions of the electron motion such as nonlinear 
oscillators p9l - l3l| or scattering off the plasma potential 
p3 [ are more appropriate. Here we do not consider such 
strong field effects. 



potential Eq. p3)) . The combination of the first two 
fields may be considered a fast harmonic perturbation 
on the latter field. Although UCPs behave entirely 
classically [H, a quantum mechanical description of 
this situation proves best suited to calculate the RF 
energy absorption. Accordingly, the electrons occupy 
bound states in the plasma potential, and can change 
states by absorption or emission of an RF photon. The 
quantum mechanical calculation of the absorption is 
given in detail in Ref. psj . In summary, perturbation 
theory is applied, in which the transition probability of 
electrons between any pair of states is given by Fermi's 
Golden rule [45|. The number of RF photons absorbed 
by the plasma equals the difference between the number 
of electron transitions to a higher state and those to a 
lower state, and the absorbed RF energy is this amount 
multiplied by the photon energy. Exploiting, in addition, 
the fact that the system dimension a is much larger 
than the typical De Broglie wavelength of the electrons, 
one can adopt the quasi-classical or Bohr-Sommerfeld 
theory to approximate quantum mechanical quantities 
by their classical analogues [1^. Although results for a 
general three-dimensional potential are available [28| . 
we will use the one-dimensional analogs because then 
the mathematics is much more transparent. This does 
not represent a major error since the energy transfer 
from the RF field to the plasma proceeds via electrons 
that move partially resonant with the applied field. This 
means that only one component of the electron trajecto- 
ries, namely the one that is parallel to the applied field, 
contributes to the RF absorption, so that the problem is 
essentially one-dimensional. Explicit calculation of the 
RF absorption in both the full three-dimensional and 
corresponding one-dimensional hard wall potential (28| 
confirms that the latter captures the general behavior. 

Expressing as before the absorbed RF power Pp due to 
coUisionless absorption in terms of an effective frequency 
I'p, it is found that (28j 



Pp IvpUp^ 



E 



\dn/de 



■ exp 



ksTe 



(14) 
(15) 



B. RF absorption by electrons in a general 
potential 

We now proceed to calculate the coUisionless RF 
energy absorption by an UCP, taking account of the 
smooth plasma potential shape shown in Eq. ()13|) rather 
than resorting to a hard wall approximation. We make 
use of the calculational method developed by Zaretsky 
et al. [1^. When forcing an UCP with an RF signal, 
the electric field in the plasma consists of the external 
RF field, the polarization field caused by any excited 
plasma modes, and the field corresponding to the plasma 



Here, r2(e) is the oscillation frequency of the classical 
trajectory a;(e, t) of a particle with energy e in the unper- 
turbed potential. 



X(e) 



2tt 



x(e,t) exp (iuit) dt (16) 



is the Fourier component of the classical trajectory at the 
frequency of the perturbation. 



Z 



exp 



de 



(17) 



is the partition function of the electron distribution over 
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the energy states, which is assumed a Boltzmann distri- 
bution here, and the sum in Eq. (jlSp is over energies 
that are roots of the equation 



{2s + l)n{es) 



(18) 



Without attempting to explain all details underlying 
Eqs. (|15m8|) here, it is noted [l^ that the only contribu- 
tions to the absorbed energy come from those electrons 
whose trajectory is in resonance with the applied field 
according to Eq. p8)) . This is the correlation between 
thermal motion and applied field also characteristic for 
bulk Landau damping. Furthermore, the contributions 
in Eq. are proportional to the spectral content 
of the trajectory at the applied frequency. However, the 
dominant frequencies in the spectrum of the trajectory 
will be on the order of the oscillation frequency ri(e). 
In a potential such as Eq. (|13l) with r replaced by x, 
this frequency will be comparable to that of a harmonic 
oscillator potential with the same curvature at a; = 0, 
that is, to 17 ~ 2Uo / {3ma'^) = wq. Therefore, it 
is expected that the RF energy absorption strongly 
depends on the ratio uj/ujq. In addition, the ratio of 
particle energy e and thermal energy kgTf, appears in 
Eq. ([T5|) . the former being limited to values smaller 
than the potential depth Uq, so there will be some weak 
secondary dependency on the ratio Uq/ {ksTf^) as well. 
These properties are indeed found below. 

In the classical UCP system the spacing between en- 
ergy levels is much smaller than the thermal energy, 
therefore the sum in Eq. (jlSp may be approximated by 
integration over s. A subsequent change of integration 
variable from s to the energy introduces an extra fac- 
tor {des/ds)~^ , which is the density of resonant states. 
This factor is obtained by differentiating Eq. (|18p with 
respect to s, yielding \dfl/de\^^^ ■ dcs/ds — 2^ juj. Ac- 
cordingly, Eq. becomes 



nmuj 





X{e) 


f 





exp 



' 2ZkBT, 
where the subscript s has been dropped 



(19) 
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FIG. 3: Model potential (red dashed line, Eq. (I20p ') compared 
to the actual UCP potential (black solid line, Eq. (|13p ). The 
parameters have been set to oji = and a such that U\ = Uq. 
The dotted parts of the model potential are not used. 



the trajectories are known analytically 
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mujix I ^ X 



(20) 



where a is a positive constant with units of length. Eq. 
((20)) is the potential of a Duffing oscillator commonly 
used to describe the motion of a mass on a cubic 
softening spring. Although this potential differs from 
the actual UCP potential Eq. we note that from a 

physical point of view the most important characteristics 
of the UCP potential are the temperature, which sets 
the potential depth J7o, and the charge density, which 
sets the curvature towq at the bottom of the potential. 
Therefore we should obtain a reasonable estimate for 
the energy absorption by choosing the model potential 
accordingly, setting the curvature muJl equal to mwg 
and the potential depth mwla^ /S = Ui equal to Uq. 
Important as well is that the infinitely differentiable 
UCP potential is modeled by an equally smooth one, and 
that both potentials approach their edge with vanishing 
slope. In Fig. [3] the two potentials are compared. 



C. RF absorption in a model plasma potential 



A particle is bound by the potential Eq. ([20)) only if 
its energy e is less than Ui. For such a bound particle 
the classical trajectory, starting at time t = at position 
a; = 0, can be shown to be given by the periodic function 

El 



Eq. (fTO)) allows explicit calculation of the absorbed 
RF power, if the classical trajectories in the potential are 
known analytically. However, for the particular potential 
Eq. (|T3|. closed expressions for the trajectories are not 
available. In order to still make quantitative estimates 
for the energy absorption, instead of Eq. (IT51) we use a 
model potential with the same general shape for which 



V 



(21) 



where sn{y, ni?) is the Jacobi elliptic function with argu- 
ment y and modulus m, and u = e/Ui is the particle en- 
ergy in units of the potential depth, and v = 1 + Vl — u. 
The frequency CI with which the particle oscillates back 
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FIG. 4: Effective collision frequency z/p for RF absorption due 
to the finite plasma size, as a function of the ratio of the fre- 
quency characterizing the potential ui to the RF frequency 
OJ. A potential depth equal to fesTe has been assumed. The 
inset shows the behavior for tJi/oi < 1 on a logarithmic scale, 
comparing the exact result Eq. (|A3|I (solid curve) to the ap- 
proximate result Eq. (|23p (dots). 



and forth in the potential is given by [4€ 



2V2K(u/w2 



-UJi, 



(22) 



where K{m?) is the complete elliptic integral of the 
first kind with modulus to. In the limit of vanishing 
particle energy e — 0, the trajectory (|21l) approaches 
harmonic motion with frequency wi, while the motion 
becomes anharmonic with the frequency monotonically 
decreasing to zero as the energy grows to Ui . 

In Appendix \X\ the absorbed power is calculated by us- 
ing Eqs. dH]) and ^ in Eq. The exact result Eq. 
(|A3[) for the effective collision frequency is plotted in Fig. 
|3]as a function of Wi/w, assuming a potential depth equal 
to ksTe- Also plotted is the asymptotic approximation, 
valid for uji/ui ^ 1, 



= C{Y)[- 

OJ J Model V 0J\ 



exp 



'V27r— 



(23) 



which fits the exact result very well. In a typical UCP, 
CT ^ 1 mm and Tg ^ 1 K [^, while in a typical RF 
experiment a;/(27r) > 1 MHz [10], so that usually the 
asymptotic regime of Eq. is in effect. The pref actor 
C (F) is given by Eq. (IA4[) and depends on the ratio 

Y = Ui/ksTe- As argued previously, the choice of model 
potential parameters that best represents the actual 
UCP potential is wi — ujq and Ui — Uq ^ ksTe , giving 

Y ^ 1. The corresponding prefactor in Eq. ([25)1 lies in 
the range C = 20 - 35 for F = 0.5 - 2.0. 



sionless RF absorption by an UCP strongly depends on 
ioi/oj, that is, on the ratio of the frequency at which the 
thermal motion of the UCP electrons takes place to the 
RF frequency. This strong dependency was anticipated 
above from the fact that the collision frequency Eq. (fTS)) 
is proportional to the spectral content of the trajectory 
at the RF frequency: when wi and lo do not differ too 
much, the RF forcing and the electron motion take place 
on more or less the same time scale, so that the electron 
motion contains an appreciable Fourier component at 
the RF frequency, resulting in resonant and efficient 
energy transfer. Since all oscillation frequencies given 
by Eq. ((22)) are in fact less than wi, the average 
oscillation frequency will be less than wi as well, so that 
in Fig. 21 the peak in the energy transfer occurs at a 
somewhat higher value than wi/w = 1, corresponding to 
a somewhat slower forcing. 

An important feature of Fig. |4] and Eq. (l23l) is the 
threshold- like behavior of Vp-. for wi > w the absorp- 
tion is significant, while for ui/lo — >■ it decreases ex- 
ponentially. The inset shows that this decrease is very 
rapid, so that collisionless absorption is completely neg- 
ligible if ciJi ^ Lo. This condition can be written as 
1 ^ uji/u] ^ oj^juj = ^2UqI (Zraa'^ufl) ~ Uth/lcw). 
Physically, this corresponds to the situation in which a 
low temperature yields by assumption a shallow potential 
with slow electrons, so that almost no electrons traverse 
the plasma within one RF oscillation. Combined with the 
lack of steep features in the smooth potential, this means 
that there is almost no electron motion available at the 
RF frequency that is susceptible to resonant absorption. 
One may thus define a critical temperature 



(24) 



From Fig. |3] and Eq. 



it is clear that the colli- 



that separates a temperature regime Tg > Tp in which 
collisionless absorption is significant and a regime Tg ^ 
Tp where it is negligibly small. Note that this behavior is 
not at all described by the hard wall approximation Eq. 
(fT2|l . The reason for this is that an electron bouncing 
between hard plasma boundaries abruptly changes its ve- 
locity at every wall collision, giving rise to high-frequency 
components essentially regardless of the velocity. There- 
fore Eq. (fT^ predicts significant collisionless absorption 
at any temperature, but is valid only for steep plasma 
potentials. 



D. Validity for the actual UCP potential 

As we just described, the collisionless absorption rate 
in the model plasma potential exponentially decreases 
with the ratio wth/(ca;). Since the physical arguments 
leading to Eq. (p4)) are valid for any general smooth 
plasma potential, also in actual UCPs the collisionless 
absorption rate will quickly decrease once the electron 
temperature is below the critical temperature Tp. How- 



ever, one may still ask whether the decay constant of this 
decrease (i.e. the factor \f2-K in Eq. (|23|) ) is also repre- 
sentative for actual UCPs, or depends on the potential 
shape. Lacking analytical expressions for the trajecto- 
ries x(t) in the UCP potential, this cannot be verified 
by explicit calculation. Nevertheless, the decay constant 
can be calculated by quantifying the asymptotic behav- 
ior of the Fourier coefficients of the trajectories, using 
the so-called Darboux's Principle [13]. This however re- 
quires considering the analytical continuation of a:(i) to 
the complex t-plane. The details are rather technical and 
are relegated to Appendix [Bl The main result is that the 
quantity |-'^(e)|^ in Eq. (fTOt for the UCP potential con- 
tains an extra factor of approximately exp (2a;/ajo) as 
compared to the case of the model potential, indepen- 
dent of the particle energy e and for sufficiently large 
Lo/ijjQ. Including this extra factor in the result Eq. (l23l) . 
the asymptotic rate of decrease of the collision frequency 
is approximately equal to 



mcr 01 = knTf 



uj J UCP 



(X exp 



^/2^ - 2 - 

/ Wo 



(25) 



Although the decay constant v^tt— 2 is smaller than that 
of Eq. (f23| and Fig. |4j it is still of the same order of mag- 
nitude. Also in the UCP case, therefore, the collisionless 
absorption is negligible for ^ w, or equivalently for 
temperatures below Tp given by Eq. p4|. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we considered three mechanisms by 
which an RF field influences the temperature of an UCP. 
First, RF energy is absorbed through the well-known 
process of collisional absorption, in which electrons gain 
energy during Coulomb collisions with ions. Second, 
the RF field modifies the TBR rate by tunnel-ionizing 
electrons from intermediate high-lying Rydberg states. 
Third, resonance between the motion of electrons in 
the plasma potential and the RF field may give rise to 
collisionless energy absorption. For all of these processes, 
naive extrapolations from well-known formulas are inad- 
equate for UCPs or strong RF fields. For example, the 
electron-ion collision frequency Eq. ^ is much smaller 
than the Spitzer frequency for strong RF fields, suppress- 
ing the collisional absorption rate. As we indicated, this 
is because the quiver velocity effectively takes over the 
role of the thermal velocity, or equivalently, because the 
temperature is replaced by the ponderomotive potential 
in the collision frequency. Likewise, the TBR rate in 
strong RF fields is much smaller than expected from 
the commonly used Te ^^^-scaling because the pondero- 
motive potential replaces the temperature in the TBR 
bottleneck energy level. Figure [5] schematically shows 
the various heating regimes in terms of the RF field 
amplitude and frequency; the strong-field effects apply 
to the area above the slanted line. As discussed in the 





(c) 



















Up= ksTe 



FIG. 5: Heating regimes for RF-driven UCPs in terms of the 
applied frequency uj and field strength Eq. (a,b): Collisional 
absorption rate according to Spitzer collision frequency and 
TBR rate according to Te ^^^-scaling; (c,d): Collisional ab- 
sorption rate according to collision frequency Eq. Q and 
TBR rate according to Eq. (a,d): Collisionless absorp- 

tion relevant; (b,c): Collisionless absorption negligible. 



previous section, collisionless absorption is only relevant 
at sufficiently high temperatures or low frequencies, as 
is represented by the area to the left of the vertical line 
in Fig. [5] 

Let us conclude by giving two numerical examples. 
The RF experiment of Fletcher et al. [l^] was well in 
the weak-field regime {a,b) of Fig. [5] according to the 
reported experimental values. Using these values in Eqs. 
(I2][8]), (|T4l) and p3| gives absorption rates per electron 
of Pci/kB = 3 K/iJS and Pp/kB = 0.002 K//^s at the 
highest reported frequency and amplitude. Considering 
the electron temperature of 100 K and the typical plasma 
expansion time of microseconds, these low absorption 
rates will not infiuence the plasma temperature and 
expansion much. For somewhat larger RF amplitudes, 
however, the collisional absorption starts to become 
significant on the time scale of the plasma expansion, 
which may be related to the high-field effects observed 
in the experiment. 

As an example in the regime (c) of Fig. ([5]), consider 
an applied field with an amplitude of 0.1 MV/m at a 
frequency of 28 GHz, which is currently available psj . 
We deliberately choose this relatively high frequency 
because otherwise the oscillation amplitude of the 
plasma electrons would exceed the plasma size at such a 
large field strength, which situation is outside the scope 
of this paper. Choosing further a = I mm, Te = 1 K 
and n = 10^ cm^^ Eqs. (Ull), (in]),(Il4l) and ^ give 
P^i/kB =4-102 K/fis and R/ue = 9 • 10"'^ fj.s-\ while 
the collisionless absorption rate is vanishingly small. 
Thus collisional absorption is expected to heat the 
plasma to the 100 K scale during the expansion time of 
the plasma, while the chance that an individual electron 
recombines is very small. Now compare these numbers 
to the corresponding results obtained from standard 
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expressions. Using the Spitzer colhsion frequency 
instead of Eq. (g]) would give Pei/ks =4-10^ K//is, 
which would predict immediate heating of the UCP to 
conventional eV plasma temperatures. According to 
the usual Te~^^^-scaling (Eq. p3| with Up replaced by 
ksTe), the TBR rate per electron would be R/ue — 50 
/is~^. Assuming an energy release of ~ kgTf. per 
recombination, this would result in a heating rate 
per electron on the order of 10^ K//is due to TBR 
alone, although of course this rate would be quickly 
quenched as the electron temperature rises. Based on 
the hard wall approximation Eq. (I12p with v — Vth, the 
coUisionless absorption rate would be Pp/ks — 1 • 10'^ 
K//is rather than exponentially small. From these 
numbers it is clear that it is essential to properly take 
into account strong field effects on the one hand, and 
the smooth UCP plasma potential on the other hand. 
For the application of a very strong microwave field to 
an UCP, it changes the predicted effect from destroying 
the plasma immediately to only heating it up moderately. 

In summary, we have analytically studied well-known 
plasma heating mechanisms and specialized them to the 
system of an UCP driven by a uniform, and possibly 
strong, RF field. Benchmarking our results against 
molecular dynamics simulations will yield valuable 
additional insights, and will also identify any additional 
RF effects that are not addressed in this paper. Among 
these are, for example, plasma cloud deformations ex- 
pected when the electron oscillation amplitude becomes 
comparable to the plasma size, relativistic effects, and 
development of plasma modes and other instabilities. 
Experiments in which RF fields are used to probe 
plasma resonances rely on adequate modeling of the 
UCP expansion dynamics, which will benefit from 
detailed knowledge of RF heating mechanisms such as 
those discussed in this paper. Furthermore, in virtue of 
comparable coupling parameters, RF-driven UCPs may 
be seen as millimetre-sized scale models of laser-driven 
solid state density plasmas. Understanding the ways 
in which ultracold plasmas interact with RF fields is 
therefore also relevant for such high-density systems. 



Substituting Eq. ^TS\\ in Eq. ([T6)) . and comparing with 
it follows that 



X{e) 



— — csch 



UJl 



K 1 



(A2) 



Using this quantity in Eq. and changing the inte- 

gration variable to u — e/C/i, results in 



I'pot 



=ttW2 — 

U! \UJl 



(A3) 



Y- 



csch^ \-^J^K (1 - u/v^)] cxp {-Yu) du 



lo ^ -'^/^ K (u/w^) cxp (— Ku) du 

where Y — Ui/ksTe. The integrations are over ener- 
gies smaller than the potential depth, corresponding to 
bound electrons, since transitions to the continuum do 
not give rise to energy increase of the ensemble that is 
left behind. When wi/w <C 1, to a good approximation 
cschZ « 2exp {—Z) in the numerator of Eq. (|A3[) . Fur- 
thermore, the argument Z of the csch-function is small- 
est at M = 1, so that the region around the upper inte- 
gration limit will give the dominant contribution to the 
integral in Eq. (jA3[) . and Z may be approximated by 
its Taylor series around u = 1. This gives cschZ w 
2 exp [-(7r/y2)(a;/wi) (1 -1-3(5/16)], where S = 1 - u. 
Similarly, in the integral in the denominator of Eq. (jA3l) , 
the elliptic function diverges at it = 1, so that again the 
region around the upper integration limit will give the 
dominant contribution, and the elliptic function may be 
approximated by its asymptotic value [i^l- This gives 
^-1/2 X {u/v"^) ~ — In ((5/64) /4. With these approxima- 
tions, the integrals in Eq. (IA3|) can be solved analytically, 
yielding Eq. ([23]), with 



C{Y) 



2567r 



y2 



3 EinF + 61n2(expr- 1) 



(A4) 



Here, Ein denotes the modified exponential integral [42|, 
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Appendix A: Effective collision frequency 

The Fourier series of the trajectory (|21[) equals p£ 



x{e,t)^2a—y2 



sin[(2n+ 1)^] 



'^i „=o sinh 



(Al) 



We use the following theorem [47|: 
The coejjicients of the Fourier series ^ a„ sin (nClt) 
of a 21T /^-periodic function y{t), which is infinitely 
many times differentiate , decay asymptotically as 
a„ cx exp(— rirn). The constant r equals min|Imij|, 
where tj denote the singularities of the function y{t) in 
the complex t-plane. 

Writing uj = (uj/fl) ■ fl in Eq. (|16|) shows that X is 
essentially the oj / fi-th Fourier coefficient of the function 
x{e, <), so that according to the theorem the integrand in 
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FIG. 6: Branch cuts (black solid lines) of the integrand of 
Eq. (|B2|) in the complex z-plane, using u — 1. In the shaded 
sectors |arg2;| < 7r/4 and |7r — argz| < 7r/4, the error func- 
tion behaves as erf (z/a) — ^ 1 as \z/a\ — >■ oo. Two possible 
contours from the origin to infinity are shown. 
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FIG. 7: Decay constant r in Eq. (|23p as a function of particle 
energy. Solid line: analytical result for the model potential 
Eq. H13p assuming cji = uq and Ui — Uo; dots: numerical 
result for the UCP potential ((20)) 



x(e, t) does not have any algebraic branch points. 



the collision frequency Eq. (jl9p is proportional to 

(X exp(-2a;r) ; t = mm\lmtj\ (Bl) 

for large oj. This expression is easily checked for 
the model potential: the elliptic function in the tra- 
jectories Eq. (PT|) has singularities along the lines 
Imt = ±u}i^ y/2/vK (l — u/v^) = ±T in the complex 
t-plane 0| . Substitution in Eq. (jBl[) yields the behavior 
of \X\^ for large w, which coincides precisely with what 
is found in Appendix [XI Eq. (jA2l) by explicit calculation. 

Applying Eq. (jBip to the actual UCP potential re- 
quires explicit expressions for the trajectories x{e,t), 
however these are not known. Instead, the inverse func- 
tion t(e, x) may be obtained by integration of the equa- 
tion of motion mcPx/dt^ = —dU{x)/dx, yielding 

pm f dz , 



Here, the initial conditions a; = and dx/dt — ^J2e/m 
at <: = have been assumed, and U {z) denotes the UCP 
potential Eq. (fT3)) with r = z. Equation (|Bip requires 
knowledge of the singularities tj of the functions a;(e,t), 
which may be categorized as either poles, logarithmic 
branch points or algebraic branch points. (More patho- 
logical singularities such as exp(l/z) at z = are not 
considered here.) An algebraic branch point in x(e,t) 
corresponds to a critical point in the inverse function 
t(e,x), at which dt/dx — 0. Differentiating Eq. (jB2p 
with respect to x, it follows that U (z) must diverge at 
such a point if the derivative dt/dx is to vanish. But 
the UCP potential Eq. (|13p is an entire function, so 
that this does not occur for any finite complex z, hence 



Considering next poles and logarithmic branch points 
in x(e,<), at such points the position diverges while the 
complex time has some finite value. In terms of the in- 
verse function Eq. (jB2p then, there exist contours Cj in 
the complex a;-plane from the origin to infinity such that 
t(e, x) tj with \tj \ < oo as a; — > oo along Cj. In view of 
Eq. (jBip we are interested in the contour that yields the 
time tj with the smallest imaginary part. A complication 
in finding this contour is the presence of the square root in 
Eq. (IB2p . because of which the integrand has branch cuts 
in the complex z-plane. Adopting the standard choice of 
letting the branch cuts coincide with the points at which 
the argument of the root is real and negative, these cuts 
start at the zeros of the function e — U{z) and extend to 
±ioo without crossing. Fig. [5]shows the resulting branch 
cut structure for the case e = Uq/2] the integrand in the 
lower half-plane is the complex conjugate of that in the 
upper half-plane. Also drawn are two possible contours 
from the origin to infinity. Now, the potential U (z) in Eq. 
(IB2p contains the error function erf (z/cr), which has the 
property that its value is close to unity for \z/a\ > 1 
in the shaded sectors in Fig. [SI while its amplitude grows 
supcrcxponentially as z — )> oo in the non-shaded sectors. 
Therefore the integrand in Eq. (jB2p will be essentially 
constant along parts of contours that cross the shaded 
sector, such as C2, so that a large contribution to the 
integral is accumulated along these parts. Hence we may 
expect that the contour yielding the smallest possible 
value of tj is the contour that avoids the shaded sectors 
altogether, that is, the contour Ci along the imaginary 
axis. With this conjecture, we calculate r in Eq. (|Bip by 
integrating Eq. (|B2p along Ci for several values of the 
particle energy e. The result is shown in Fig. [3 together 
with the analogous result for the model potential. As is 
clear from the figure, for any particle energy r for the 
UCP potential is approximately one unit w'^^ less than 
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that for the model potential. Hence, asymptotically for model potential, independent of e. The resulting rate of 
large w, the quantity in Eq. will contain an decrease of the collision frequency is given in Eq. (1^51) . 

extra factor exp {2uj /ujq) as compared to the case of the 
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